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Abstract Bound states of complex scalar fields (boson stars) have long been proposed 
as possible candidates for the dark matter in the universe. Considerable work has al- 
ready been done to study various aspects of boson stars. In the present work, assuming 
a particular anisotropic matter distribution, we solve the Einstein-Klein-Gordon equa- 
tions with a cosmological constant to obtain bosonic configurations by treating the 
problem geometrically. The results are then applied to problems covering a wide range 
of masses and radii of the boson stars and the relevant self interaction parameters are 
calculated. We compare our results with earlier treatments to show the applicability 
of the geometrical approach. 

Keywords Exact solution • Einstein-Klein- Gordon equation • Boson star • Dark 
matter. 



1 Introduction 

According to recent cosmological observations, approximately 96% of the total material 
content of the universe is exotic in nature out of which 73% is gravitationally repulsive 
(dark energy) and the remaining 23% is attractive in nature but exists in the form 
of dark matter T,'2l,'3l. The quest for such exotic matter in the universe remains an 
outstanding problem in cosmology and astrophysics. In some models, massive axion 
of mass m ~ 10^^ eV is considered as a standard dark matter particle (with no self- 
interaction). These massive axions may collapse to form compact objects known as 
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axion stars of masses comparable to the mass of a planet. The survival of such axion 
stars has been questioned by Seidel and Suen[4]. Moreover, these small objects cannot 
explain the observed abundance of dwarf galaxies and the dark matter density profiles 
in the core of galaxies. The view that dark matter in the galactic halos could be due 
to the presence of scalar field confogurations has been considered by many[5ll6], (see 
also Ref.[7] for a recent review). Hu et aZ[8] showed that the large scale structure of 
the galactic halo could be explained by dark matter composed of ultralight scalar 
particles of mass m ~ 10~^^ eV. Lee and Koh[9] have investigated boson stars with a 
self-interacting scalar field as the model for galactic halos. These results are consistent 
with the observation that the mass of a scalar particle should be greater than 10~^^ eV 
not to disturb radial stability 10_. 

The large abundance of dark energy in the universe is also a field of active research. 
It appears that the dark energy could be well described by a cosmological constant [llj. 
In this paper we will accept this interpretation. The effect of the cosmological constant 
on the mass or radius of a compact boson star will be negligible. However, we will work 
in the general framework where the cosmological constant will be retained [12] • 

Historically, the concept of a gravitationally bound state of scalar particles, was 
proposed by Kaup[T3] and Ruffini and Bonazzola[14]. Since a boson star is prevented 
from gravitational collapse by Heisenburg uncertainty principle, the maximum mass of 
a boson star with no self-interaction, is very small. However, Colpi et ttZ[T5] showed that 
a self-interaction of the scalar particles could yield compact objects of masses compa- 
rable to neutron stars. In the absence of direct observational evidence of fundamental 
scalar particles constituting a boson star, various models have been proposed to study 
the gross features of such hypothetical objects P^17lfT51ll9ll2UII21ll22) . One of the main 
motivations to study these objects is to understand its possible role in explaining the 
dark matter in the universe, although boson stars made up of Higg's particles, dilatons 
and other scalar particles may also exist. The choice of the potential corresponding to 
the self-interaction is very important in the calculations of such models. A quartic type 
potential {A<1>'^) is, in general, used for the self-interaction [T5]|21) . though some works 
on higher order self-interactions [22] and different forms such as cosh-type potential[l] 
I23j have also been reported in the literature. The motivation for choosing a certain 
potential form is to obtain configurations which are consistent with astrophysical and 
cosmological observations. For example, assuming a cosh-type potential of the form 
V{<P) — Vo[cosh(A\/87rG#) — 1], Alcubierre et al i23^ showed that one can construct 
a self-gravitating scalar object or oscillaton of critical mass M^rit ~ lO"'^'^ Mq which 
is comparable with the dark matter content of a galactic halo where, the parameters 
have values A ~ 20.28, Vb ~ (3 x 10"^"^ mp/)"*, ~ 1.1 x 10"^^ eV. Henriques and 
Mendes[TS] considered a combined boson-fermion star and showed that it is possible 
to generate a wide variety of configurations ranging from objects of atomic sizes and 
masses of the order of 10~^^ Mq to objects having galactic masses of the order of 
10^^ Mq and radii extending upto a few light years. Of particular interest is a con- 
figuration where the fermion core and the bosonic halo could form a supermassive 
compact object of mass M ~ 10"'^^ Mq for the bosonic particle mass mf, ~ 10""^^ GeV 
and fermionic particle mass mf ^ 10~^ GeV. The method that one usually follows 
in such calculations is that one makes a choice for the potential and then solves the 
Einstein-Klein-Gordon system numerically by fixing the values of the coupling constant 
and the value of the scalar field at the origin. It will be interesting if we assume the 
mass and radius of a boson star first and then look for the potential corresponding to 
the assumed configuration. In the present work, we explore this alternative approach. 
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where the mass of the scalar particle and the coupling constant get determined by 
the field equations. Essentially we will solve the Einstein-Klein-Gordon system with 
an energy momentum tensor corresponding to a self-interacting complex scalar field in 
the presence of a cosmological constant. Current observational studies done by High- 
Z Supernova Team[2] and Supernova Cosmological Project group[3] clearly indicate 
that, though small, we nevertheless need a positive cosmological constant to explain 
the current acceleration of the universe expansion. The current value of the cosmolog- 
ical constant is estimated to be ylc ~ 3 x 10~^^ cm~^[24]. In our model, it helps us 
to construct a star with a proper boundary. In general, in the standard calculation of 
a boson star, one imposes the condition that the scalar field vanishes asymptotically 
which makes it difficult to define the radius properly (usually defined as where 90% 
of the mass is contained) . However, a cosmological constant term makes it possible to 
define the radius in a conventional way (as in a fermion star), i.e., where the radial 
pressure vanishes. Similar observations may also be found in some other papers [3ll25]. 
In a recent paper^G] the radius of the dark halo is described as the distance where pres- 
sure is equal to the cosmological constant Ac of the order of magnitude 10~^^ gm/cm^ 
or alternatively, the cosmic microwave background radiation (CMBR) which is about 
10""^^ gm/cm'^ and the corresponding radii are 1 Mpc for Ac and 100 Mpc for the 
CMBR. The exterior of the halo is described by the Schwarzschild de-Sitter spacetime. 
Note that boson stars have already been studied in de-Sitter as well as anti-de-Sitter 
universe [27]. 

Our paper is organized as follows. The boson star model is developed in Section 2 
where, by making use of an ansatz for one of the metric potentials for an anisotropic 
matter distribution, we generate a solution for the Einstein-Klein-Gordon system. In 
Section 3, we utilize this solution to obtain different boson star configurations and 
determine the nature of the self interaction. We discuss our results in Section 4 and 
outline the physical implications of our studies along with some comments on the earlier 
work on the possibility of the formation of a boson star[28] and on its stabilitv[29] under 
radial perturbations. 



2 Boson star configuration 

We write the energy-momentum tensor for a scalar field in the form 

T^J = i + " l9^J [g"^'' $*m<l^n + Vm^)] , (1) 

where ${r,t) is a complex scalar field and V^(|^|^) is the potential of self-interaction of 
the scalar field. The scalar field is assumed to be of the form 

#(r,t) = 0(r)e~'"', (2) 

which guarantees a spherically symmetric static matter distribution. We write the 
metric for the spherically symmetric space-time in the standard coordinates 



(3) 



where 7(r) and ^{r) are to be determined. The Einstein's field equations (with c — 1) 
are then obtained as 
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where vie is a cosmological constant term. A prime (') here denotes differentiation with 
respect to the radial coordinate r. The Klein-Gordon equation 



for the line element ((3)l takes the form 
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Eq. (|4])-((6]) together with Eq. ((8} comprise the system of equations describing the boson 
star. 

Note that by analogy with an anisotropic fiuid, we may identify Eq. Q-© as 
density(p), radial pressure(pr) and tangential pressure(px) equations, respectively [10]. 
It then follows that in the interior of a boson star, pressure is anisotropic due to the 
term (f)''^e~'^^ . Clearly pr > pj_ in this model. If Zi = 8nG{pr — pj_) is the measure of 
pressure anisotropy, i.e., 



then Eq. ([5]) and Eq. ((6} may be combined to yield 
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which is a non linear second order differential equation with two coupled functions 
7(r) and /i(r-). To solve Eq. (I10|l we first make use of an ansatz for one of the metric 
potentials given by Vaidya and Tikekar 30 in the form 
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where R and k are dimensionless parameters. Eq. (|10p then gets the form 

(1 + fe - fcx^)^,, + kx9^ + k{k + 1)^ + \/ ' ^ = 0, 

(1 — x^j 
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where we introduced the following transformations 
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To solve Eq. (|12p the radial dependence of the anisotropic parameter A needs to be 
specified. We make a choice for the anisotropic parameter[31j 



where a > is a constant. Eq. (I12|) then gets the form 



(1 - z^)^zz + zl'z + [fc(l + a) + 1]^ = 0, 



(13) 



(14) 



where z = yjk/{k + l)x. Using the properties of Gegenbauer function and Tschebyscheff 
polynomial, a general solution of Eq. (|14|l may be obtained as [32) 
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where, /3 = \Jk{l + a) + 2, C, = cos~^ z, and A and 5 are constants. The energy density 
(p), radial pressure {pr), tangential pressure (px) and the anisotropic parameter [A) 
are then obtained as 
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Thus all the physical parameters can be obtained once the geometry of the star is 
specified. 

Combining equations. ([9]), (|lip and (|19p and integrating we also obtain a solution 
for the scalar wave function in the form 



a = (70 + 
where ctq is an integration constant. 



/ ak . _i 
— sm z. 
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2.1 Determination of the Potential 



To analyze the form of the potential, we now express Eq. (|4])-(|6| and Eq. ([8} in dimen- 
sionless forms as 
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where we have made the following rescaling: r* = rm, 

Lu/m, yl* = ylc/m^, p = 8nGp, pr = STvGpr and pj^ = SvrGpx- Here m is the mass of 
the scalar particle and a prime now denotes differentiation with respect to r*. 
We rewrite the Eq. (f24)) in the form 
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where n = (fc + l)/k. Defining the potential in dimensionless form as 
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we write Eq. (|25|) in the form 
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where, i7 = QRif. The physical quantities are then obtained as 
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Combining the two sets of equations, i.e., Eq. ([T6)l - ([T9)l and Eq. lf28|l - IpTjl . we 
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where ylc = 9>-kGR^ Ac- Eq. (|32|l - (|35| ) thus help us to obtain our physical quantities 
from two different perspectives - one from the grometry and the other from the energy- 
momentum part of the field equations. 
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2.2 Boundary conditions 

To calculate the physical quantities of the geometrical part of the field equations, we 
impose the following boundary conditions. 

(1) At the boundary of the star (r — b) the interior solution is matched to the external 
Schwarzschild-de Sitter metric which gives 

e^-(-^) - (l - ^ - ^) , (36) 



We note that the effect of Ac for highly compact objects like neutron stars is neg- 
ligible in this model. However, for large objects whose compactness is very low, the 
contribution due to Ac can not be ignored. 

(2) The radial pressure pr should vanish at the boundary which gives 

where = {k/{k + 1))(1 - b'^/R^). 

For given values of mass and radius, Eq. (|36p - p8p may be used to calculate the 
values of the constants A, R and 5 for assumed values of the curvature parameter k 
and the anisotropic parameter a. Thus all the physical quantities can be obtained once 
the geometry is specified. 



3 Numerical results 



Using the results of Section 2, we now choose different bosonic configurations and 
calculate the interaction parameters. The method is the following. 

For a given mass M and radius 6, we calculate the constants from the boundary 
conditions, viz., R (from Eq. S (from Eq. ^) and A (from Eq. We 

combine Eq. (|32p and (|33|l to obtain 
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where there are two unknown parameters, fi and ap. We plot (p+pr) against the radial 
coordinate r from the centre to the boundary of the star using the geomerty part (left 
hand side) of Eq. (|39}. By suitably choosing the values of H and (jq, we fit this curve 
with the one obtained from the left hand side of Eq. (|39p . This has been shown in 
Figure [T] 

We now choose the scalar wave potential in the form 
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a \ 5 \ A \ -an \ Q \ R* \ A \ Vq \ Vo{Mev/im^) \ m(oV) 


Case I: fe = 100, M = IMq, b = 15 km and R = 305 km 


0.5 


2.738 


40.094 


6.993 


29.710 


24.463 


1.599 


42.2 


27.4 


1.6 X 10-^^ 


1.0 


2.644 


48.478 


9.965 


24.948 


18.168 


0.648 


60.8 


39.5 


1.2 X 10"^^ 


2.0 


2.409 


62.550 


14.229 


23.179 


9.710 


1.572 


90.7 


58.9 


6.3 X 10-^^ 


Case II: k = 100, M = IMq, b = 10 km and R = 155.68 km 


0.2 


2.647 


31.133 


4.462 


33.320 


28.546 


6.352 


27.8 


69.3 


3.6 X 10"^^ 


0.5 


2.615 


35.962 


6.943 


23.651 


19.670 


2.083 


38.1 


94.8 


2.5 X 10"^^ 


1.0 


2.511 


43.231 


9.899 


19.733 


14.090 


1.055 


54.8 


136.5 


1.8 X 10"^^ 


1.5 


2.386 


49.661 


12.196 


18.536 


9.786 


1.392 


69.6 


173.5 


1.2 X 10"^^ 


2.0 


2.257 


55.446 


14.149 


18.234 


3.638 


11.305 


81.6 


203.4 


4.6 X 10-^^ 


Case III: k = 100, M = 1.5Mq, b = 10 km and R = 113.25 km 


0.5 


2.378 


29.789 


6.822 


14.787 


10.679 


4.348 


34.6 


162.9 


1.8 X 10"^^ 


1.0 


2.252 


35.337 


9.738 


12.079 


5.874 


4.083 


49.1 


231.7 


1.0 X 10"^^ 


1.2 


2.195 


37.384 


10.700 


11.592 


3.185 


10.789 


54.8 


258.0 


5.6 X 10"^^ 


Case IV: fc = 2, M = IMq, b = 10 km and R = 28.58 km 


0.5 


2.301 


1.099 


0.207 


3.590 


2.006 


7.12 


1.66 


122.73 


1.4 X 10-^^ 


1.0 


2.274 


1.243 


0.376 


2.928 


1.740 


2.49 


1.94 


143.43 


1.2 X 10"^^ 


2.0 


2.157 


1.511 


0.693 


2.621 


1.089 


2.17 


2.58 


190.75 


7.5 X 10"^^ 


2.2 


2.128 


1.561 


0.753 


2.604 


0.907 


2.90 


2.70 


199.62 


6.3 X 10-^^ 


2.5 


2.085 


1.635 


0.840 


2.590 


0.529 


8.11 


2.87 


212.19 


3.6 X 10"^^ 



Table 1 Results of different bosonic configurations. 



where the constant potential term Vq, the coupling constant A and the mass of the scalar 
particle m are yet to be fixed. To this end, we rewrite the potential V in dimensionless 
form as 

V = 4ttGR'^V ^Vo + rU^ +ARla'^, (40) 
where A = XjAixGrv? and Vb = /I-kGR^Vq. Differentiating Eq. (gOl) we get 

'^^IRia + ^ARiG^. (41) 
da 

We now plot ^ against a using Eq. (f27)l and Eq. (|4T|l as shown in Figure [2l This 
fixes the values of 7?* and A. The mass of the scalar particle is then obtained from the 
relation m = Rt, j R. 

Finally making a suitable choice for the constant Vq, the energy density p and the 
radial pressure pr are plotted against the radius r as shown in Figure [3] and Figure 
21 respectively. The variation of radial pressure with density (EOS) both from the 
geometric part and the matter part are shown in Figure [S] 

Thus we get a complete description of the boson star. Following the technique, 
we have considered a wide variety of bosonic configurations and the results have been 
compiled in Table [1] 



4 Discussion 

The geometric approach developed in this paper helps us to study the form of the 
potential of a wide variety of bosonic configurations. The key features of our studies 
are given below: 
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2 4 6 8 10 

r 

Fig. 1 (p+pr) plotted against radius r (km). The solid curve originates from the geometry and 
the dotted curve is derived from the matter part of the field equations. We took, M = 1 Mq, 
fe = 10 km and a = 1. 




0.38 0.4 0.420.440.460.48 0.5 



a 

Fig. 2 plotted against a. The solid curve originates from the geometry and the dotted 
curve is derived from the matter part of the field equations. We took, M = 1 Mq, b = 10 km 
and a = 1. 

— Considering a boson star of mass M = 1 Mq and radius 6 = 10 km, we have shown 
the variation of scalar field a against radius r for three different values of the 
anisotropic parameter a as shown in Figure [6] The scalar field decreases smoothly 
towards the boundary from the centre and is nodeless, as desired. Possible radial 
excitations of the scalar field configurations have not been considered here. 

— For given values of M, b and k, the mass of the scalar particle m decreases and 
the constant potential term Vq increases with the increase of the anisotropic factor 
a. The coupling constant A first decreases, remains almost constant for a specific 
range of a, and then increases. These variations are shown in Figure [T] 

— Note that a cannot be made arbitrarily large. In our model, the maximum value 
of a decreases with increasing compactness for a given k. However, for the same 
compactness it increases if k is made small. 
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Fig. 3 Energy density p plotted against radius r. The solid curve originates from the geometry 
and the dotted curve is derived from the matter part of the field equations. We took, M = 
1 Mq, fe = 10 km and a = 1. 




r (km) 

Fig. 4 Radial pressure pr plotted against radius r. The solid curve originates from the ge- 
ometry and the dotted curve is derived from the matter part of the field equations. We took, 
M = 1 Mq, 6 = 10 km and a = l. 



— In this model we have considered a quartic type potential only, however, higher 
order terms in the potential may also be considered. Variations of the potential V 
against the scalar field a and with respect to r inside the stellar interior have been 
shown in Figure [8] and Figure [Q] respectively. 

— Our method can be used to study a wide variety of bosonic configurations. In Table 
[2] we have compiled some of our results which are comparable with the results 
obtained by earlier workers. 

— Boson star of mass of the order of a planet: The model can be used to 
describe a boson star of mass of the order of a planet. For example, if we consider 
a star of mass M = 10~^ Mq and radius b = 10"^ km, the mass of the scalar 
particle turns out to be m ~ 10~^ eV. In CDM model, it is observed that massive 
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Fig. 5 Radial pressure pr plotted against density p. The solid curve originates from the 
geometry and the dotted curve is derived from the matter part of the field equations. We took, 
M = 1 Mq , 6 = 10 km and « = 1. 

axions of similar order of masses can collapse to form compact object of mass 
M ~ 0.6mp2/m ~ 10"^ M© g]. 

— Axion star: Stable axion stars have masses less than 0.846 Mq and radii greater 
than 20.5 km with mass of the scalar field m ~ 10"^*^ eV, close to the lower bound 
of axions [33]. For a star of mass 0.8 M© and radius 20.5 km, we get the mass of 
the scalar particle in the same order. 

— Boson star in galactic centre: Some galactic centres are supposed to contain 
super-massive compact dark objects. It is believed that these ultra-compact objects 
are either supermassive black holes or very compact clusters of stellar size black 
holes [34]. Observational data seem to allow the possibility that these could well be 
bosonic stars. In a boson star, the radius b > M {m~pf)^E\. Data collected for the 
central object of a galaxy by Genzel et al'^ show that for a boson star model as 
in Ref. [35], the mass and radius of the central object should be M ~ lO^M© and 
b ~ 10^ km, respectively. The mass of the corresponding scalar particle turns out 
to be m ~ 10~^^ eV. These results are consistent with our calculations. 

— Boson star of galactic scale: In this model, we find that for a scalar particles 
of mass m ~ 10 eV and constant potential Vq ~ 10 Mev/fm , a boson 
star of mass M ~ 10^^ M© and radius b ~ 10^^ km may be obtained. These 
values are comparable to the scalar field dark matter model of Ref. [23], where 

~ 1.1 X 10"^^ eV and Vq ~ 10"^^ Mev/fm^ for a potential of the form 
V(<?) = Vo[cosh(AV87rG'?) - 1]. The size of the galaxy of mass ~ 10^^ M© was 
about 10^^ km. Thus the results are in good agrement with the present model. 

— We have shown that a wide variety of boson stars may be obtained in this model 
by either scaling or by considering different compactness. Obviously not all of these 
stars will be realistic. One also needs to consider the stabifity of such configurations 
under radial perturbations. This may provide constraints on the mass of the scalar 
field[9], and therefore, the mass of the boson star. For radial stability, the mass 
of boson should be m > 10~^* eV[9lll0j. In this model, with such scalar particle 
masses, one can describe a boson star of mass and radius as high as A4 ~ 10^* M© 
and b ~ 10^^ km, respectively. 



12 



M (Mq) 


fc(km) 


Vo (Mev/fm^) 


A 


m (eV) 


10-" 


10--'' 


9.619 X 10" 


2.08 


2.46 X lO-'' 


0.8 


20.5 


8.4 


1.34 


10-11 


10^ 


10^ 


9.619 X 10-11 


2.08 


2.46 X 10-1''' 


1012 


1013 


9.619 X 10-23 


2.08 


2.46 X 10-23 


1018 


1019 


9.619 X 10-35 


2.08 


2.46 X 10-29 



Table 2 Variations of scalar particle masses for different boson star configurations with a = 
0.5 and k = 100. 




Fig. 6 Scalar field a plotted against radius r (km). Dotted, solid and dashed curves are for 
a equal to 0.5, 1 and 1.5 respectively. 
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a 

Fig. 7 Variations of A (soUd Une), Vq/IO in Mev/fm^ (dotted Une) and m/lO-n in eV (long 
dashed line) with a. 
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120 




0.35 0.4 0.45 0.5 
a 

Fig. 8 Scalar field potential V plotted against scalar field a. Dotted, solid and dashed curves 
are for a equal to 0.5, 1 and 1.5 respectively. 



300 
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Fig. 9 Scalar field potential V plotted against radius r. Dotted, solid and dashed curves are 
for a equal to 0.5, 1 and 1.5 respectively. 

To conclude, although the present model is based on a particular geometry, depend- 
ing on two parameters (Vaidya-Tikekar model), and a simple anisotropic distribution, 
the resulting configurations, at least a class of them, may be physically relevant. This 
is confirmed by comparing the results of our geometrical approach with earlier results, 
obtained by direct numerical integration. Whether boson stars can actually account 
for the dark matter, or at least a part of it, will require extensive analysis of the 
observational results of the CMB radiation, the Supernova data and other relevant 
astronomical results. The stability of the boson star configurations under radial per- 
turbations also remains an interesting issue. For scalar fields without self-interaction, 
extensive numerical calculations done by Seidel and Suen |28) have shown interesting 
stability properties of the configurations for finite perturbations. Whether boson stars 
with self-interacting fields have similar behaviour remains to be checked. It may be also 
be pointed out that Seidel and Suen[29] have also suggested a mechanism, the gravita- 
tional cooling, which permits the scalar field to get rid of the excess kinetic energy so 
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that a bound configuration can be formed (instead of a virialised cloud). These results 
make the possibility of the formation of boson stars look more likely. We hope to take 
up these issues elsewhere. 
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